Dendritic flux avalanches and nonlocal electrodynamics in thin superconducting films 
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We present numerical and analytical studies of coupled nonlinear Maxwell and thermal diffusion 
equations which describe nonisothermal dendritic flux penetration in superconducting films. We 
show that spontaneous branching of propagating flux filaments occurs due to nonlocal magnetic 
flux diffusion and positive feedback between flux motion and Joule heat generation. The branching 
is triggered by a thermomagnetic edge instability which causes stratification of the critical state. 
The resulting distribution of magnetic microavalanches depends on a spatial distribution of defects. 
Our results are in good agreement with experiments performed on Nb films. 

PACS numbers: 74.20.De, 74.25. Qt, 74.25.Fy 



Penetration of magnetic flux in a type-II superconduc- 
tor can result in nonequilibrium pattern formation, such 
as magnetic macroturbulence |l|, kinetic front rough- 
ening [3, magnetic microavalanches and dendritic 
structures Dendritic flux penetration has been re- 
vealed by magneto-optical imaging (MOI) on multiple 
scales *~ 1 — lOO^m much greater than intervortex spac- 
ings in YBa 2 Cu30 7 |3, Nb IflEzI, Nb 3 Sn||, and MgB 2 0. 
Similarity of these dynamic flux patterns in different ma- 
terials indicates a generic collective behavior of vortices. 

Recently it has been shown both experimentally and 
theoretically that dendritic flux penetration is due to 
a positive feedback between moving flux and the Joule 
heating coupled b y a highly nonlinear voltage-current 
characteristic |j, |[J, . The resulting thermal bistability 
of current-carrying superconductors gives rise to switch- 
ing waves between a cold superconducting phase and a 
hot resistive phase self-sustained by the Joule heating 
Dendritic flux penetration in superconductors has 
analogs in the theory of pattern formation out of equilib- 
rium |l2j and instability of solidification fronts yjj- 

Dendritic flux dynamics in superconductors was ob- 
served in numerical analysis of nonisothermal magnetic 
diffusion in a slab in a parallel field, when flux penetra- 
tion was triggered by a local heat pulse 0] . However, ex- 
periments have been mostly done on films in a perpendic- 
ular ramping magnetic field in which case magnetic flux 
diffusion becomes st rong ly nonlocal due to long-range in- 
teraction of vortices 14]. In this Letter we calculate den- 
dritic flux penetration controlled by nonlocal magnetic 
flux diffusion coupled to thermal diffusion in thin films. 
We report a novel nonlocal mechanism of flux branching, 
which captures salient features of dendritic flux penetra- 
tion in superconducting films. 

We consider a thin film strip of the width w along the 
y-axis and thickness d <C w in the xy plane perpendicular 



to the magnetic field Hq. Distributions of the magnetic 
induction, B(r,i), and temperature T(r,t) are described 
by the Maxwell equation coupled to the heat diffusion: 

Cd t T = VkVT- {T-T )h/d + JE(J,T), (1) 
8 t B = -VxE(J,T), VxH = J(S(z). (2) 

Here C(T) is the heat capacity, k(T) is the thermal con- 
ductivity, h(T,To) is the heat transfer coefficient to the 
coolant or substrate held at the temperature To, and 
E = JE(J, T)/J is the electric field, which strongly de- 
pends on T(r, f) and the sheet current density J(r, t). 

The E(J,T,B) characteristic accounts for a resistive 
flux flow state with E = (J — J c )pf for J > J c and a 
low-resistive flux creep state with E = E c cxp( J — J c ) / J\ 
for J < J c , where J C (T, B) is the critical current density. 
We use the following interpolation formula expressed in 
terms of observable parameters ficj : 



E = p F J\ Ml + exp( J - J c )/Ji], 



(3) 



where J\(T) is logarithmic flux creep rate ( J\ -C J c below 
the irreversibility field B < B*), and Pf{T) = p n B/B c2 
is the flux flow resistivity. 

We consider weak Joule heating, for which the most es- 
sential temperature dependence comes from E(T), while 
other parameters may be taken at T = Tq. The rela- 
tion between current and the z-component B z in a film 
is given by the non-local Biot-Savart law. Expressing 
J x = d y g and J y = —d x g in Eq. © in terms of the cur- 
rent stream function g(x,y,t), we obtain the equations 
for g and the dimensionless temperature 9: 



rg = K[d x [r(3, W*9 + d y [r(j, 6)d y g] - rH (t)] , (4) 

e + afr(j,V)- (5) 



V 2 6 



Here we define the operator K in the Fourier space, 
K = ^ k sin(k x x) shi(k v y)gk/k, where k x = irn/L and 



mn/w with integer m and n to ensure zero nor- 
mal component of J at edg es of a rectangular film of 

Furthermore, 8 = (T — 
3 = J/Ji 



width w and length L 

T )/(T* - T ), J C (T*) = 



eclg 



H z {x,y) = J2k ex P( — dk — ikr)g^k/2 + H (t), the factor 
exp(— dk) accounts for a finite film thickness, the deriva- 
tives in Eqs. (@J and (JSJ are taken over normalized time 
t/th and coordinates r/Lh, and the nonlinear resistiv- 
ity r(j,6) = ln[l + exp(j — j c (0))]/j is obtained from 
Eq. ©. Here the thermal length Lh = {dn/h) 1 / 2 and 
time th = Cd/h define the spatial scale and the cool 
down time of T(r, t) at a frozen J(r). Hereafter we take 
W) = io(l - h for 61 < !) and ic = for 6» > 1, and 
Ji (T) = const , assume that a uniform magnetic field 
H (t) = H t is ramped up with the rate H$, and neglect 
the field dependence of J c . Evolution of 6*(r, t) and g(r, t) 
is controlled by two dimensionless parameters: 



p.J 2 d 



2p F C 



h(T* - T ) ' 



(6) 



Here r = t m /th is the ratio of magnetic and thermal 
diffusion times, and a quantifies the Joule dissipation. 
Magnetic nonlocality strongly reduces r = Tod/2Lh in a 
film as compared to tq = pon/ppC in the bulk. Indeed, 
flux diffusion over a distance Lh along a film takes t m ~ 
dLh/D m , while thermal diffusion takes £/, ~ L\jDh, 
where D m — (Xo/pF and = n/C. For Nb films with 
d = 0.5pm, at 4.2K, (k ~ 0.2W/cmK, ft ~ lW/cm 2 K, 
Pn ~ 3 • lCT^cm 0, and C ~ 2 • 10- 3 J/cm 3 K), we ob- 
tain L h = {dn/h) 1 / 2 ~ 0.03mm, t h = Cd/ft ~ I0- ? s, 
(i/L/j ~ 10 -2 , t ~ 0.1, with r decreasing as T increases. 

We used Eqs. (QJ, © to study flux penetration in 
a film with periodic boundary conditions (b.c.) along 
x and d y 8 = g = at y = Q,w. Eqs. J2J|, were 
solved numerically by a quasi-spectral method based on 
the Fast Fourier Transform; up to 1024 x 512 harmonics 
were used. To implement non-periodic b.c. in the y- 
direction we used the domain of doubled length with the 
condition g(x,y) = —g(x,2w — y). Calculated steady- 
state distributions of B z (x) are very close to those of 
the Bean model |l7). In the majority of numerical runs 
for t <C 1, we observed spontaneous avalanches induced 
by ramping magnetic field Hq(€) = Hot, starting from 
a zero field cooled state. We also took into account 
randomly-distributed macroscopic defects modelled by 
je(r) = j [l - -J2i ftcosh^dr - r;|/£ )] where q { , r, 
and £o determine the strength, the position and the ra- 
dius of the i-th defect. Extended current-blocking defects 
have been commonly revealed by the MOI [lfij |. 

We start with flux penetration in a film with no macro- 
defects, Fig.^and Movie 1 in [T^|. The ramping magnetic 
field first caused penetration of a stable cold flux front. 
Then an instability causing periodic modulations of tem- 
perature and propagation of hot magnetic filaments over 
the preceding smooth magnetic flux distribution devel- 
ops at the film edge. Once the first wave of magnetic 




FIG. 1: Temperature distributions in a film with no defects 
for H = 5,h/t h , t = 0.0025, a = 0.008, system size 600L h x 
1501/h (only quarter part is shown) and jo = 20 at t = 16.75th 
(a) and t = 17.75th (b); (c) magnetic flux pattern for t = 20th 



filaments reaches the central line where the magnetiza- 
tion currents change direction, the filaments widen and 
start splitting at the ends, similar to that of a slab in 
a parallel field ^(|- At the same time, a second wave of 
hot filaments start propagating from the film edge mostly 
between the paths of the filaments of the first wave. The 
new filaments are wider and exhibit shape distortions due 
to interaction with preceding filaments. Eventually the 
film cools down and the flux penetration stops, resulting 
in a frozen multi-filamentary structure |l7j . 

To address the mechanism of the edge instability, we 
performed a linear stability analysis of steady-state dis- 
tributions T(y) and Eo(y). Because E(J,T) depends on 
J and T, linearization of Eqs. J3J and J5J gives two cou- 
pled equations for the perturbations 59 and Sg, which 
depend on Eq(jj) and T{y). The problem can be simpli- 
fied by the fact that the flux stratification period i in Fig. 
1 is much smaller than the film width w. Because both 
Eo(y) and T(y) vary on the scale ~ w, we can neglect 
slow variations of Eo(y) and T(y), taking uniform T = Tq 
and Eq — Bqw/2 at the film edge. Then the solutions, 



d y 5T 



at 



which satisfy the boundary conditions SJ y 
the film edge (y = 0), take the form 66 cx e xt+lkx cos qy, 
and Sg cx e xt+tkx s\x\qy. From Eqs. Q, (J5J, we obtain 
the following dispersion relation: 



(A + 1 + q 2 + k 2 



-(3)(\T l ^k 2 + q 2 + q 2 + sk 2 ) 

+<z 2 /?(i 



T 

«) = 0, 



(7) 



where A and (k,q) are measured in t h , and L h , re- 
spectively. The dissipation control parameter (3 = 
(dJ/h)dE/dT depends on E a : if E(J) = E c exp[{J - 
J c )/Ji] then (3 ~ EoJcdldJc/dTl/^h for J < J c . The 
parameter Tj = /ioV dhn/2pC ' , is similar to that in Eq. 



3 



|jBJ, except that pp is now replaced with the differential 
resistivity p(E) = OE/dJ, and s = E/Jp(E) ~ Ji/J c is 
the flux creep rate. 

Eq. J7J describes coupled thermal and magnetic diffu- 
sion modes, the factor \Jk 2 + q 2 accounting for the mag- 
netic nonlocality. Positive eigenvalues A(fc, q) correspond 
to unstable modes resulting in spontaneous thermomag- 
netic structures 0|- F° r slow flux diffusion r ^> 1, ther- 
mal perturbations with A~/3 — f — fc 2 — q 2 are unstable 
above the thermal runaway threshold, j3 > 1. In this case 
A is maximum at fc = q = so no periodic structures are 
expected. A completely different situation occurs for fast 
flux diffusion r <C I for which Eq. J7J yields 



A = /3-l 



2 9 2 /3(l + a) 



fc 2 - 



sfc 2 



(8) 



The spectrum of q is determined by the full set of bound- 
ary conditions in the theory of flux jumps For 
further qualitative analysis we take q ~ tt/26, where 
6J x (x,b) = 0, b = [I - I/cosh(B /-B p )W2 is the 
width of the flux penetrated critical state region, and 
Bp — p J c /ir |20|. For a given q, the increment X(q,k) 
passes through a maximum at the wave vector k m , which 
defines the period I — 2n/k m of the fastest growing 
thermomagnetic structure along the film. Here sfc 2 ^ = 
[q 2 /3s(l + s)] 1 / 2 — q 2 , thus X(k,q) is maximum at the fi- 
nite k m if /3 > /3j = q 2 /s(l + s), or B > Bi. For s < I, 
we obtain B; and ^ in normal units: 



B, = 



16b 2 



Ab 3 \dJ c /dT\ : 



(B a /B t y/ 2 - 1 



(9) 



The period I decreases as 1/s and J3n increase, in agree- 
ment with our numerical results |17| . The branching 
instability with k m > and Re A(/3 c ,fc m ,g) > oc- 
curs at /? = (3 C in Eq. (jHJ), that is, B > B c = 
h,hP c /dbJ c \dJ c /dT\, where 



fc' 2 = (l + 2?) 1 ' 2 +qy/TTT/i- 



(10) 



For r,; < 1, Eq. (JTUI) gives /3 C > ft for all = n/2b(t), 
so a thermomagnetic structure with the finite period 
1{(3 C ) "C 6 develops as the width b{t) exceeds b c , where 
Bq = B c (b c ). Notice that Eq. (JJJ) defines a region 
t\ < t < T2 in which A is complex, which manifests 
itself in temporal oscillations of growing flux structures. 

Next we consider dendritic flux penetration initiated 
by macroscopic defects, both at the film edge and in the 
bulk. Such defects can trigger local flux jumps even if the 
critical state in the bulk is stable j2]| and cause branch- 
ing instability of flux filaments in a slab in a parallel field 
|l0| . Selected results of flux patterns in a film arc shown 
in Fig. 2. For j a = 1/s = 20 in Fig. 2a, edge defects 
produce flux fingers superimposed with a smooth flux 
front. This behavior is characteristic of any superconduc- 
tor with a highly nonlinear E(J), for which a defect of 



size £o produces a much larger disturbance ~ £o/s across 
the current flow [2^] . The flux fingers widen and split at 
the ends as they collide with the central line where mag- 
netization currents change direction [lOj . As jo increases, 
hot flux filaments in Figs. 2b and 2c get thinner and 
start branching even before they reach the central line. 
Then new fingers start growing between the defects due 
to the edge instability considered above, see animations 
for more details 0|- Moreover, in wider samples (Fig. 2 
c, d) the finger undergoes multiple branching giving rise 
to the characteristic flux dendrites. This new branching 
mechanism, which was not observed in simulations of flux 
patterns in a slab [lfj , is principally due to nonlocal flux 
diffusion in films. Thinning the filaments as s decreases 
follows from Eq. © , while the branching shape instabil- 
ity is facilitated by magnetic nonlocality, which results in 
a weaker damping of short-wave electromagnetic modes 
A m cx fc as compared to local flux diffusion for which 
A m cx fc 2 . Another new effect at higher values of jo is a 
"giant" flux avalanche in Fig. 2c, which starts propagat- 
ing from the region with no surface defects after the first 
wave of smaller flux filaments reached the center. 

Fig. 2d shows flux patterns in a film with randomly- 
distributed bulk defects. In addition to the branching due 
to magnetic nonlocality, propagating flux filaments can 
undergo splitting caused by local transient heat spikes as 
they collide with defects. This results in local shape in- 
stability of the filaments and their subsequent branching 
similar to that obtained for local magnetic flux diffusion 
flf)j |. With the decrease of jo (which is equivalent to an 
increase of To or Hq in the experiment) flux filaments 
become wider and eventually start overlapping, forming 
a continuous flux front. However, even in this case a sig- 
nificant front roughness still persists both due to a size 
distribution of individual filaments and due to local heat 
releases as the flux front collides with defects 17]. 

The dendritic flux penetration can be regarded as 
avalanches of vortex bundles which do not trigger a global 
flux jump in the whole sample. Such avalanches pro- 
duce local temperature spikes, partial fluxpenetration 
and a step on magnetization curves M(B) Nonlo- 
cal electrodynamics of films strongly facilitates branching 
flux propagation, which requires faster magnetic diffu- 
sion t = tq y/ dh/n <C I characteristic of thin films with 
dh k. Notice that because b ~ w(Bq/ B p ) 2 for small 
Bo, both critical ramp rates Bi cx Bq and B c cx Bq 6 
strongly depend on the applied field B for B Q <C B p but 
level off for Bq > B p as b — > w/2. Furthermore, both Bi 
and B c increase as To increases. Thus, for a given Bq, the 
branching occurs at lower T above a certain field Bi(T), 
in agreement with many experiments 0, 0, S @ • 

Fig. 3 shows MOI of flux penetration in two different 
Nb films in ramping fields. One 9mmxl.8mmx0.5fim 
Nb film described in detail elsewhere exhibits flux pat- 
terns similar to those in Figs. 2a, b where dendritic flux 
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FIG. 2: Flux penetration in a film of 320L h x 40L h with 30 
random edge defects with £o = 0.25Lh, and amplitudes qi 
uniformly distributed between and 0.4, for Hq — 3Jj/th, 
a = 0.08, r = 0.0025 and t = 25t h for j = 20 (a) and 
jo — 80. Black and white correspond to the Meissner and 
vortex phases, respectively; (c) Flux penetration in wide film 
of 600L h x 150L h for H = 5Ji/t h , j = 60, 20 edge defects. 
The "giant" avalanche develops in a defect free region; (d) 
Flux pattern in a film with 500 randomly distributed bulk 
defects at jo = 80. See also Movies 2-6 in |17j . 




FIG. 3: MOI of flux branching in Nb films of Ref. 
at T=6.2K and B =31.2mT (a) and 4.7K and B =36.5mT 
(b). "Giant flux avalanche" at 4.5K in Nb film of @ (c). See 
also Movies 7-8 in [ill ]. 

penetration is initiated by surface defects. For this film 
(k = 0.2W/cmK and h — lW/cm 2 K), we obtain dh/n ~ 
10 -4 and the thermal length L h = (dn/h) 1 / 2 = 0.03mm, 
much smaller than the film width w — 1.8mm (i.e. the 
film width is about 60Lh as in Figs. 2a, b), in which case 
the magnetic nonlocality does play the key role. Fig. 3c 
shows MOI of an effectively wider (4mmx4mmx0.1/xm) 
film [2^ (w ~ l0 2 Lh), in which in addition to small 
microavalanches near the film edge a "giant" avalanche, 
similar to that triggered by a laser pulse in YBa2Cu307 
H and those in Fig. 2c, d develops. 

In conclusion, we proposed a mechanism of flux frag- 
mentation in superconducting films caused by coupling 



of nonlocal flux diffusion with local thermal diffusion. 
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